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I. 


INTRODUCTION 


This report summarizes results obtained through the support of NASA Grant 
NAGW-2412. The objective of this project is to conduct a comprehensive investigation of the 
radio wave emission observed by the Planetary Radio Astronomy (PRA) instrument on board 
Voyager 2 as if flew by Neptune. This study has included data analysis, theoretical and numerical 
calculations, ray tracing, and modeling to determine the possible source mechanism(s) and 
locations of the Neptune radio emissions. We have completed four papers, which are included 
in the appendix. The paper "Modeling of whistler ray paths in the magnetosphere of Neptune" 
investigated the propagation and dispersion of lighting-generated whistler in the magnetosphere 
of Neptune by using three dimensional ray tracing. The two papers "Numerical simulations of 
bursty radio emissions from planetary magnetospheres" and "Numerical simulations of bursty 
planetary radio emissions" employed numerical simulations to investigate an alternate source 
mechanism of bursty radio emissions in addition to the cyclotron maser instability. We have also 
studied the possible generation of Z and whistler mode waves by the temperature anisotropic 
beam instability and the result was published in "Electron cyclotron wave generation by 
relativistic electrons". 

Besides the aforementioned studies, we have also collaborated with Drs. Sawyer, King, 
Romig, and Warwick of the PRA team to investigate various aspects of the radio wave data. 
Two papers have been submitted for publication and the abstracts of these papers are also listed 
in the appendix. 

In the next section we summarize the research results. 

II. SUMMARY OF RESEARCH RESULTS 

A. Modeling of Whistler Ray Paths in the Magnetosphere of Neptune 

During the Voyager 2 flyby of Neptune, 16 "whistlerlike" events near the magnetic 
equator were observed. Using a simple empirical plasma model based on observational constraints 
and an offset titled dipole magnetic field model, a three dimensional ray tracing study is 
performed to investigate the propagation and dispersion of lightning-generated whistlers in the 
magnetospheres of Neptune. Source positions for the whistlers are chosen along magnetic field 
lines with L < 3 in agreement with observations. Whistlers with frequencies of about 10 kHz 
generally propagate very nearly the resonance cone angle for most of the ray path, which 
produces dispersions considerably larger than those predicted for quasiparallel propagation. In 
general, the whistlers were found to propagate about 10 degree in magnetic longitude during one 
hop and to cross L shells. These results strongly support the existence of whistlers at Neptune. 
The large dispersions are attributed to the large cold plasma population at Neptune, which leads 
to propagation very near the resonance cone angle for a large part of the ray path. 


B. Numerical Simulations of Bursty Radio Emissions from Planetary 
Magnetospheres 

The Voyager 2 spacecraft has observed both smooth and bursty radio emissions 
from Uranus and Neptune. These emissions are known to be freely propagating primarily in the 
right-hand circularly polarized (RCP) mode with the bursty emissions having duration periods 
as short as a few tenths of a second and the smooth emissions being observed over periods of 
a few hours. While the smooth emission is probably due to the cyclotron maser instability, some 
other processes might be responsible for the generation of the bursty emissions. It is proposed 
that one important difference in mechanisms is that the smooth emissions are associated with 
continuous injection of electrons while the bursty emissions are associated with impulsive 
injection. In the latter case, the electron distribution can develop a beam feature with a 
temperature anisotropy. It is shown via one-dimensional (three velocity) relativistic particle 
simulations that while the beam may be initially unstable to an electrostatic instability, this 
instability quickly saturates and eventually the beam is unstable to a strong electromagnetic beam 
instability which utilizes the temperature anisotropy as free energy; a beam feature is not 
explicitly needed for the growth of this electromagnetic instability. The radiation generated by 
the instability is able to convert particle energy to wave energy at similar levels at the maser 
instability. It is also postulated that some of the radiation generated below the local X mode 
cutoff may also be able to escape and be detected remotely via mode conversion between regions 
where field-aligned currents produce local perturbations in the magnetic field. 

C. Electron Cyclotron Wave Generation by Relativistic Electrons 

Besides the freely propagating extraordinary (X) and ordinary (O) mode waves, Z mode 
waves with frequencies below the electron cyclotron frequency have also been observed in 
planetary magnetospheres. It is shown that an energetic electron distribution which has a 
temperature anisotropy, or which is gyrating about a DC magnetic field, can generate electron 
cyclotron waves with frequencies below the electron cyclotron frequency. Relativistic effects are 
included in the stability analysis and are shown to be quantitatively important. The basic idea 
of the mechanism is the coupling of the beam mode to slow waves. The unstable electron 
cyclotron waves are predominately electromagnetic and right-hand polarized. For a low density 
plasma in which the electron plasma frequency is less than the electron cyclotron frequency, the 
excited waves can have frequencies above or below the electron plasma frequency. This 
instability might account for observed Z mode waves in the magnetosphere of Neptune. 
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A three-dimensional ray tracing program is used to investigate the propagation and dispersion of lightning-generated 
whistlers in the magnetosphere of Neptune. A simple empirical magnetosphenc plasma model based on observational 
constraints is used for the plasma density distribution, and an offset tilted dipole model is used for the magnetic field. 
Source positions for the whistlers are chosen along magnetic field lines with L < 3 in agreement with observations. 
Whistlers with frequencies of about 10 kHz generally propagate very near the resonance cone angle for most of the ray path, 
which produces dispersions considerably larger than those predicted for quasiparallel propagation. Dispersions computed 
using our "standard*’ plasma density model arc over 30,000 s VHz for "one-hop" whistlers which is comparable to observed 
values. In general, the whistlers were found to propagate about 10° in magnetic longitude during one hop and to cross L 
shells. These results strongly support the existence of whistlers at Neptune. The large dispersions are attributed to the large 
cold plasma population at Neptune which leads to propagation very near the resonance cone angle for a large part of the 
ray path. 


Introduction 

During the Voyager 2 flyby of Neptune, which occurred on 
August, 25, 1989, the plasma wave instrument detected 16 
"whistlerlike" events near the magnetic equator at radial distances 
ranging from 1.30 to 1.99 R N [Gurnett et al „ 1989, 1990]. The 
frequency ranged from 6.1 to 12.0 kHz, and the dispersions fit the 
Eckersley law for lightning-generated whistlers. The existence of 
whistler mode emissions (which propagate at frequencies less than 
the minimum of either the gyrofrequency, f# or the plasma 
frequency, f p ) sets a limit on the density at the point of 
observation. The magnetic field was directly measured at the 
point of observation, and the local gyrofiequency during the 
observations was at least 100 kHz. Belcher et al. [1989] have 
reported that plasma density measurements near closest approach 
(CA) are typically 0.01 to 0.1 cm' 3 with a maximum of 1.4 cm 3 
near the magnetic equator. These values are dependent on ion 
composition which is uncertain. However, the instrument was 
probably not sensitive to cold plasma because the Faraday cups 
were not favorably oriented at the time of CA (R. McNutt, private 
communication, 1990). We also note that Barbosa et al. [1990] 
have estimated the density at this time based on plasma wave 
measurements of ion cyclotron harmonics to be in the range of 1 
to 10 cm \ again dependent on the ion composition. Sawyer et 
al. [1990] reported what appear to be electron cyclotron 
harmonics at this time in the Planetary Radio Astronomy data, 
which may also indicate higher densities based on the cyclotron- 
maser instability. The topic of plasma density was discussed at 
great length by Gurnett et al. [1990], and led the authors to 
suggest that the signals were in fact whistlers. Gurnett et al. 


believe the density at the point of observation must be at least 30 
to 100 cm* 3 . 

One result that is difficult to understand is the extremely large 
dispersions, typically 26,000 s VHz. The dispersion is defined by 
the Eckersley law, t = DHf where t is the travel time, / is the 
frequency and D is a contant called the dispersion. For low 
frequencies and propagation parallel to the magnetic field the 
dispersion constant D is given by 

D = <l/2c) (1) 

where ds is elemental distance along the ray path. Based on 
minimum and maximum estimates of the density in the 
magnetosphere, Gurnett et al. calculated that path lengths ranging 
from -50 to over 2000 R** are required to achieve the observed 
dispersion. While whistlers observed at Earth are in rare cases 
known to have bounced as many as 200 times, it is difficult to 
imagine an electromagnetic wave surviving such long path lengths 
without suffering significant damping, scattering, or other 
processes that would strongly degrade the amplitude. 

In this paper we present a preliminary study of whistler 
propagation based on three-dimensional ray tracing. We 
emphasize that these are initial results and more work needs to be 
done to understand the details of the observations. The results 
show that very large dispersions are obtained even without 
multiple bounces. The large dispersions are the result of 
propagation very near the resonance cone angle for a large part 
of the ray path. Such dispersions are considerably larger than 
predicted by the Eckersley law, which is not valid at large wave 
normal angles. 
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Plasma and Magnetic Field Models 

The magnetic field model is an offset tilted dipole which is a 
refinement (OTD2) of an earlier model developed by Ness et al. 
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[1989]. This model is based on a dipole, with magnetic moment 
of p = 0.13 GR*„ tilted by about 45° with respect to the rotation 
axis and is offset from the center by 0.55 R N . This model is 
known to be quite accurate along die spacecraft trajectory at 
distances greater than about 3 but may have significant 
deviations from the real field close to the planet (N.F. Ness, 
private communication, 1990). The large offset produces a strong 
asymmetry in the field strength between northern and southern 
magnetic hemispheres, and the tilt produces a strong longitudinal 
asymmetry. As will be shown, these asymmetries have noticeable 
effects on whistler propagation. 

The plasma model consists of an ionospheric model and a 
magnetospheric model. The ionospheric model is due to 
Shinagawa and Waite [1989]. Specifically we have used the 
model represented by curve "a” in their Figure 1, page 946. This 
was their "standard case" model and produced peak densities of 
about 6 x 10 4 cm' 3 at an altitude of about 1700 km. This model 
differs from the results of the Voyager radio science observations 
[Tyler et al 1989]. The maximum density obtained from the 
radio occultadons was about 2.5 x 10 3 cm' 3 at an altitude of about 
1400 km, but it is to be noted that these observations were made 
near the terminator. Most likely the electron densities are much 
higher on the dayside of Neptune, particularly near die subsolar 
point. We used a simple linear fit to the logarithmic function of 
Shingawa and Waite as follows: 

= 10W "‘ (2) 

with m = -3.16 x 10 4 , n t = 13.18, and r the radial distance in 
kilometers. For the magnetospheric electron density we have 
introduced a simple empirical model with a radial dependence 
only: 

n = n, + (3) 

with C = 99.91, k = 0.567, n a = 0.09, and R N = 24,765 km. The 
constants were selected to produce a density of about 30 cm' 3 at 
L = 3. This value is approximately midrange between the 
estimate of the density range suggested by Gumett et al . [1990]. 
For this initial study no local time dependence was incorporated 



DISTANCE ALONG HELD LINE (R n ) 

Fig. 1. Density versus distance along a magnetic field line. The footprint 
of the magnetic field line chosen was X = 40°, A = 140°. The curves are 
for equation (4) with A = 0 and A = 1. 


in die density model. We have also examined a modification of 
this magnetospheric model by introducing a latitudinal 
dependence with a peak density at the geographic equator. The 
specific form of the latitudinal dependence is given by 

n' = n[ 1 + Acos(A)] (*) 

where n' is the density, A is a variable, and X is the geographic 
latitude. A plot of the density versus distance along a sample 
magnetic field line for both equations (3) and (4) is shown in 
Figure 1. 

Ray Tracing 


Ray tracing calculations were performed using the three- 
dimensional code that has been discussed in the past [cf. Menietti 
et al „ 1990]. This code is based on the cold plasma theory of 
Sttx [1962] and Haselgrove’s [1955] set of first-order differential 
equations. Early in the study it was discovered, however, that 
whistler wave propagation near the resonance cone was occurring 
because of the high ionospheric and magnetospheric densities. 
When this happens, the index of refraction becomes so large that 
numerical instabilities prohibit further integration. To circumvent 
this problem, a modification of the basic equations was introduced 
on the basis of the scheme of Cerisier [1970], also utilized by 
Cairo and Lefeuvre [1986]. The fundamental change was the 
introduction of the following logarithmic derivative: 


_£jn(cosfy-cosfy ) * 
ds 


2sin\\tcos\\f^L - jLcosfy, / (cosfy-costy f ) 
ds ds r 


(5) 


with 

Jlcosy = 1 dS - l 

ds ' (S - P) ds (5 - Pf ds ds 


where ¥ is the wave normal angle, s is the ray path, and ¥ r is the 
resonance cone angle, which according to Stix [1962] is given by 

tan 2 * v, ) = -P/S ( 6 ) 


with P and 5 defined by Stix [1962], As discussed in detail by 
Cerisier [1970], because of the slow logarithmic dependence, 
equation (5) allows the evaluation of ¥ and numerical integration 
near the resonance cone angle. 

The Eckersley law dispersion given by equation (1) is only 
valid for low frequencies and small wave normal angles, such that 
flf e « cos(y). This assumption is called the quasi-parallel 
approximation [Helliwell, 1965]. The quasi-parallel 
approximation fails for propagation near the resonance cone 
angle, particularly when the ray path is at larger distances (r ~ 
2 R N ) from Neptune. As we describe below, we initially calculate 
the time delay, T, for a large number of rays at / = 10 kHz in 
order to obtain candidate source positions. To compute T we use 
the formula T = / dslv g where v g is the group velocity. We later 
estimate D by comparing T for a number of ray paths, from the 
same source point, at different frequencies. The initial wave 
normal angles are adjusted such that the ray paths for each 
frequency intersect one point in the magnetosphere. 


Results 

The whistlers observed by Voyager 2 all occurred within about 
35° of the magnetic equator and at 1.4 < L < 2. The L values 
were obtained from N.F. Ness (personal communication, 1991). 
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Using the OTD2 model we identified source positions near die 
surface of Neptune for which L < 3. Generally speaking there 
were two such positions along each longitude meridian of 
Neptune, one in each hemisphere. Initially we confined our 
attention to source positions located nearer die north (weaker) 
magnetic pole because the lower hybrid frequency, was much 
lower in this case. This is true because 

WL = i ifl + i/£ + W*Q (7) 

and is strongly dependent on the local field strength through f e . 
In this equation, f d and fa are the ion cyclotron frequency and ion 
plasma frequency, respectively. As an example, for a source near 
the north magnetic pole (X = 5°, A = 140°), f d - 8.5 x 10' 3 kHz, 
£ ~ 54 kHz, f c ~ 132 kHz, and/jj, - 3.1 kHz. In contrast for a 
source near the much stronger south magnetic pole (X = -30°, A 
= 250°), ~ 3.1 x Iff 2 kHz, /„. - 54 kHz, f c ~ 1710 kHz, and f w 

~ 32.1 kHz. When / < f L H , the index -of-reff action surface is 
closed, and reflection occurs when /approaches or if the wave 

normal angle approaches 90°. A wave propagating away from the 
surface in a region where » /(where magnetic field strength 
is large) can be trapped in the ionosphere, and a wave from the 
magnetosphere can be reflected back. For/> the index-of- 
refraction surface changes to an open surface with a discontinuity 
at the resonance cone angle. 

The rays were launched at an altitude of 1700 km, near the 
peak ionospheric density, at a wave normal angle that 
corresponded to a wave vector directed radially outward (vertical) 
from the planet A vertical initial wave normal direction is used 
because whistler mode waves from lightning are expected to be 
refracted to near vertical due to the large index of refraction in 
the ionosphere. Generally, the initial wave normal angle was 
about 25° to 30°. The frequency was 10 kHz, and source 
positions were chosen every 10° of longitude. Source positions 
with the smallest local magnetic field strength usually produced 
the largest time delays. This follows approximately from the 
equation for O which shows that the dispersion is inversely 
proportional to Vf c (see equation (1)). In Table 1 we list a 
summary of results for various source positions and for two 
density models as identified in the last column. The angle 0 is 
the azimuthal angle about the magnetic field; i.e., this angle 
determined the plane in which the ray was launched. The 
approximate L shell for a particular ray path as it crossed the 
magnetic equator is designated as L lvg . This L value is typically 
different from the L value of the source point because of 
deviations of the ray path from a fixed L shell. 

The time delays listed in the table are the total integrated 
values at a point near the magnetic equator before the ray 
propagates to the opposite hemisphere. Many of the travel times 
in Table 1 are much larger than comparable values calculated at 
Earth or at Jupiter. The reason for the large travel times is 
twofold. First, the enormous dimensions of Neptune produce a 
ray path much longer than any obtained at Earth, which causes a 
corresponding increase in the dispersion. Second, the wave 
normal angles of the rays are very near the resonance cone angle 
for a large pan of the ray path. At Jupiter the magnetospheric 
densities are apparently much lower than those at Neptune, and 
the magnetic field strength is much larger. Consequently, the 
lower hybrid frequency will be larger than the whistler frequency 
for much of the ray path thereby resulting in an open index-of- 
refraction surface with no resonance cone. 

It is clear from the table that source points near a longitude of 
130° produce the largest values of T. For 0 - 58°, X = 25°, A = 


TABLE 1. Selected Source Positions and Tune Delays (f= 10 kHz) 
Source Position Density 


Model 


(J. 

deg 

K 

deg 

A, 

deg 


r. 

> 

A 

50 

-10 

70 

2.7 

11.4 

0 

49 

0 

100 

2.7 

7.6 

0 

57 

18 

130 

2.8 

3.9 

0 

67 

35 

150 

2.7 

3.6 

0 

174 

70 

260 

26 

3.8 

0 

34 

10 

10 

2.4 

6.4 

1.0 

72 

-20 

30 

1.8 

1.6 

1.0 

49 

0 

1 00 

2.7 

15.5 

1.0 

53 

10 

120 

2.9 

9.4 

1.0 

58 

25 

130 

28 

61.6 

1.0 

61 

25 

140 

28 

35.1 

1.0 

67 

35 

150 

27 

233 

1.0 

93 

60 

180 

2.7 

4.6 

1.0 

103 

65 

190 

26 

4.9 

1.0 

123 

70 

210 

24 

83 

1.0 

174 

70 

260 

25 

120 

1.0 

85 

49 

330 

22 

6.0 

1.0 


130° we obtained T = 61.6 s. This is about 1/4 of the value of T 
- 250 s reported by Gurnett et al . [1990], In Figure 2 we show 
a three-dimensional plot of the ray path for this last case. 
Superimposed on the figure is the magnetic field line along which 
the source was located. The solid circle indicates the south 
magnetic pole. Clearly the ray path diverges substantially from 
the magnetic field line and extends to large L values. Two 



Fig. 2. A three- dimensio nal plot of the ray path of the whistler for a case 
where T - 61.6 s. The field line along which the source point was located 
is shown as well as the south magnetic pole (black dot). 
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perspective plots of the ray path in two dimensions are shown in 
Figure 3. The total change in magnetic longitude is about 10°. 
In Figure 4 we plot 4* and 4 / r versus ray path, s , , and in Figure 5 
we plot die index of refraction, N, and vjc (group velocity/speed 
of light) versus s. We may use these figures to understand a little 
more clearly the behavior of the whistlers in the ionosphere and 
magnetosphere of Neptune. 

In Figure 4 we see that die wave normal angle, 4*. is initially 
about 23° and that the wave propagates along the density 
gradient The resonance cone angle is near 90° at this poinL As 
the wave propagates out of the ionosphere, the wave normal angle 
initially decreases as the wave leaves the ionosphere, and then 
increases, rapidly approaching the resonance cone angle. In 
Figure 5 we see that the index of refraction initially decreases as 
the wave propagates out of the ionosphere but soon begins to 




Fig. 3. Two perspectives in two dimensions of the same ray path shown 
in Figure 2. The top panel is the perspective from the north geographic 
pole. 


increase rapidly as 4* approaches 4 / r As propagation proceeds 
into the magnetosphere and both f p and f c decrease, 4* approaches 
4 / r and the magnitude of 4 / gradually increases and then decreases 
rapidly when 4* is no longer close to 4^ The wave again begins 
to approach Neptune and a new field line footprint in the opposite 
hemisphere. During this time both 4* and N steadily decrease. 
As die wave enters the ionosphere for the second time near the 
conjugate footprint, at the south magnetic pole, 4*, increases to - 
90°, and 4 / begins to increase. Because the south magnetic pole 
surface field strength is much larger than the north polar field 
strength, the gyrofrequency becomes relatively quite large, and 
fplfc becomes less than 1 (Figure 6). The index of refraction 
reaches a minimum and begins to increase as fjjf e begins to 
increase, once again becoming greater than 1 well inside the 
ionosphere. 



Fig. 4. The wave normal angle, 4', and the resonance cone ingle, 4 , „ 
versos ray path for the case with T = 61.6 s. 



Fig. 5. The index of refraction and vjc versos ray path for the case with 
T = 61.6 s. 
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Fig. 6. Gyrofrequency and plasma frequency versus ray path for the case 
with T = 61.6 s. 


Calculation of Dispersion 

In order to estimate the dispersion of the rays we have 
evaluated the time delay for several frequencies for specific 
source points. The technique is tedious because, for a fixed 
source point in the ionosphere we require that the ray for each 
chosen frequency intercept a given observation point in the 
magnetopshere. The observation point was selected to be near the 
magnetic equator at a distance of - 2.4 R N . This point was not an 
actual spacecraft position, but was chosen to be near sources 
which produced the largest dispersion. There are at least four 
parameters that could be adjusted at the source point to 
accomplish this: X F, p, X, and A. In the limited time and finances 
available for this study, we were able to find ray paths at three 
frequencies that intercepted at the common point for whistlers 
with a source point at X = 25°, A = 130°. We see from Table 1 
that the total time delay for this source point at a frequency of 10 
kHz is about 62 s. 

To minimize the time required to obtain satisfactory results we 
limited the variable parameters to and p. We found 
interceptions for three frequencies with at most a 1° change in p. 
The initial value of the wave vector k varied by a maximum of 
0.1° from the vertical. We display plots of 1/// versus T in 
Figure 7. The dispersion for this case is found from a least 
squares fit to be D ~ 34,800 s /Hz. The latter value is 
comparable to the value of - 26,000 s /Hz obtained by Gurnett 
et al. [1990] for most of the whistlers. It is important to note that 
Gurnett et al. [1990] obtained D from the observations by 
measuring the slope of the curve of 1 />/f versus delay time. Their 
values are strictly determined by arrival time measurements as a 
function of frequency and are not dependent on the validity of 
equation (1). Our calculated values of D are obtained by 
essentially the same technique that was used to obtain the 
observed values. 

Summary and Discussion 

We have shown that whistler propagation in the magnetosphere 
of Neptune leads to dispersions comparable to those observed by 


DISPERSION CURVE 



DELAY TIME, T (sec) 

Fig. 7. 1/7/ versus delay time T for whixtlen with a source at X = 25°, A 
= 130°. 

Voyager. The large dispersions are the result of propagation over 
large distances with wave normal angles very near the resonance 
cone angle. In the absence of damping, the waves initially 
propagate approximately along magnetic field lines and drift in 
magnetic longitude by about 10° per hop. The actual dispersion 
is clearly dependent on the density along the ray path, which is 
not known with certainty at this point For a magneto spheric 
density model consistent with that suggested by Gurnett et al. 
[1990] we have obtained dispersions of - 35,000 s /Hz for one- 
hop whistlers. This value of D is comparable to those observed 
by Gurnett et al. [1990] by the plasma wave instrument on board 
Voyager 2. Our results not only strongly suggest that the 
signatures observed by Gurnett et al. [1990] are in fact whistlers, 
but make the ’’requirement" of tens of multiple hops suggested by 
those authors unnecessary. However, there are still some 
discrepancies between our results and the observations of Gurnett 
et al. which we discuss below. 

The observations indicate that all of the observed whistlers 
have approximately the same dispersion. Gurnett et al. [1990] 
suggested that this may be due to longitudinal drift of the ray 
path of the whistlers from a localized source. We suggest as an 
alternative to this explanation that only single-hop whistlers from 
a localized source may exist Our results indicate that most of the 
dispersion has occurred by the time the wave has traveled to the 
magnetic equator. Thus a satellite observing whistlers from a 
localized source would see a relatively constant dispersion. This 
does not preclude the existence of multiple-hop whistlers, 
however. 

Another factor that we must consider is the damping of the 
waves. We can obtain a crude estimate of the possibility of 
Landau damping by equating the local ionospheric thermal 
velocity to the parallel phase velocity of the whistler. We do not 
expect Landau damping outside the ionosphere because the 
plasma temperature would be much lower. For an ionospheric 
temperature of about T = 950° [Tyler et al 1989] we find the 
critical index of refraction for the onset of Landau damping is 
given by 

Merit - 2.5 x 10 3 cos( \|/ ) (8) 
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In the ionosphere of the opposite hemisphere, the ray of Figures 
2 and 3 yields ¥ ~ 42°, which yields N ^ > 1500 while the 
calculated value is N ~ 25. For all points along the ray path the 
largest wave normal angle of the wave is ¥ - 63° (Figure 4), 
which yields > 1100. The largest index of refraction of the 
wave for all points along die ray path is N - 90. We conclude 
therefore that no Landau damping is expected 

Most of the whistlers reported by Gumett et al . [1990] were 
observed near the magnetic equator at a longitude of about 250° 
with L < 2. Our calculations indicate that the largest dispersions 
should result from whistlers with sources near geographic 
coordinates of 20° < X < 40° and 130° < A < 150°. For these 
sources the ray path intercepts a longitude of -170° near the 
magnetic equator, and with 25 < L tvg < 3.0. Attempts to make 
the model whistlers intercept the region of observation yielded a 
time delay of less than 10 s for one-hop whistlers. 

From Figure 4 and Table 1 of Gurnett et al. [1990] the total 
travel time for the observed whistlers in the frequency range of 
9 < /< 11 kHz is seen to be between about 230 s and over 300 
s. For all of the cases of this study the total travel time of the 
waves was less than 100 s, even though the model dispersions are 
comparable to those observed This would indicate that the actual 
ray path or plasma parameters have not been exactly reproduced 

It is our belief that while in this initial study we have not been 
able to reproduce the observations exactly with regard to location 
we have shown that extremely large dispersions are to be 
expected for whistlers at Neptune. These large dispersions are 
primarily due to a large cold plasma density resulting in 
propagation near the resonance cone angle over great distances. 
By modifying the surface magnetic field model and/or the 
magnetospheric and ionospheric density model, it is quite possible 
to produce model whistlers that agree with the observations with 
respect to both location and dispersion. Such a study is most 
likely tedious, but may also provide an indirect measure of the 
surface magnetic field strength and cold plasma density, which at 
this point are still uncertain. 

It is also not known at this time exactly what effect higher- 
order moments of the magnetic field might have on whistler 
dispersion. We are interested in determining what parameters 
control if or how rapidly Hf approaches x ¥ r At this time we know 
that higher magnetospheric density in combination with lower 
magnetic field strength increases the time delay. The initial angle 
of the wave with respect to the density gradient, the ionospheric 
and magnetospheric density, and the curvature and strength of the 
magnetic field all control how rapidly 4* approaches The 
relative importance of each is not well known at this time. We 
have found, however, that the ionospheric density appears to be 
the least sensitive of the parameters. We thus believe that 
magnetic field strength near the source and magnetospheric 
density are the main factors in producing large delay times and 
consequently large dispersions. 
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The Voyager spacecraft observed both smooth and bursty radio emissions from Uranus and 
Neptune. These emissions are known to be freely propagating primarily in the right-hand 
circularly polarized (RCP) mode with the bursty emissions having burst periods as short as 
a few tenths of a second and the smooth emissions being observed over periods of a few 
hours. While the smooth emission is probably due to the electron cyclotron maser instability, 
some other processes must be at work to produce the bursty emissions. It is proposed that 
one important difference in mechanisms is that the smooth emissions are associated with 
continuous injection of electrons while the bursty emissions are associated with impulsive 
injection. In the latter case, the electron distribution can develop a beam feature with a 
temperature anisotropy. It is shown via one-dimensional (three velocity) relativistic particle 
simulations that while the beam may be initially unstable to an electrostatic instability, this 
instability quickly saturates and eventually the beam is unstable to a strong electromagnetic 
beam instability which utilizes the temperature anisotropy as free energy; a beam feature 
is not explicitly needed for the growth of this electromagnetic instability. The radiation 
generated by the instability is able to convert particle energy to wave energy at similar 
levels as the maser instability. For u> p */n, £ 1, most of the wave energy is LCP and is 
trapped. However, for w p ./0, £ 1 the dominant mode is RCP. The induced waves are in the 
form of modified whistlers when the beam speed is low and ~ 1 but are in the freely 

propagating x mode branch when w p «/ 0, ~ 0.4. The amount of radiation with frequencies 
above the local x mode cutoff increases with beam speed. It is also postulated that some 
of the radiation generated below the local x mode cutoff may also be able to escape the 
plasma and be detected remotely via mode conversion between regions where field-aligned 
currents produce local perturbations in the magnetic field. 


1. Introduction 

Observations of planetary radio emissions from Uranus 
and Neptune obtained by the Voyager 2 spacecraft 
have shown clear evidence of two distinct classes of 
radiation which are characterized as smooth and bursty. 
At Uranus, the smooth emission is seen as broadband 
continuous emissions, from a few kilohertz to about 
900 kHz and continuous in time over periods of a few 
hours. Many studies have suggested that the b-smooth 
emission is generated by the cyclotron maser instability 
[e.g., Kaiser et al., 1987; Menietti et a/., 1990]. The 
Uranian broadband bursty (b-bursty) emission is seen 
in the range 300-700 kHz and occurs at times which 
are intermediate between the observations of the smooth 
emissions (see Plate 1 of Farrell and Calvert [1989] for 
an excellent graphical display). Thus the smooth and 
bursty emissions are quite distinct in frequency range, 
duration, and spacecraft event time. Uranian narrowband 
bursty (n-bursty) emissions show further frequency and 
temporal fine structure [Farrell et a/., 1990]. 
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Many early studies based on straight-line propagation 
[cf. Leblanc et a/., 1987; Farrell and Calvert , 1989] 
indicated the source of the b-buxsty emission was 
colocated with the smooth emission near the south 
magnetic pole. These early studies generally agree that 
the emission propagates at large wave normal angles 
and suggest the cyclotron maser instability (CMI) as the 
source mechanism. However, Curran et aL [1990] have 
shown, via a ray tracing study, that refractive effects are 
important and that the source of the bursty emissions 
appears to be at lower latitudes for emission initially 
propagating at large wave normal angles. The conclusion 
of this latter study was that since the source region 
is different from that of the b-smooth emission, it is 
possible that the emission may be due to a distinct 
plasma distribution at lower latitudes and/or a different 
source. 

Observations of the radio emission obtained by the 
Voyager 2 flyby of Neptune also show clear examples of 
both a bursty and a smooth component [e.g., Warwick 
et al., 1989; Farrell et a/., 1990]. The bursty emission 
is observed between 500 kHz and 1.3 MHz while the 
smooth emission is seen in the frequency range f < 500 
kHz. The bursty emission is also generally observed 
at distinct times from the smooth emission, and their 
source regions appear to be different [Farrell et al ., 
1990; Ladreiter et a/., 1991]. The emission appears to 
occur in episodes with strong bursts of duration < 6 s, 
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and the polarization is consistent with the right-hand 
extraordinary mode. Generally the bandwidth of the 
Neptune bursty emission is much shorter than that of 
the b-bursty emission at Uranus, with the former having 
a bandwidth of less than 20 kHz. 

It has been suggested that both the smooth and bursty 
emissions are due to the electron cyclotron emissions, 
similar to the proposal for the emissions at Uranus [e.g., 
Farrell et al , 1990 and references therein]. However, the 
above differences in the characteristics of the smooth 
and bursty emissions suggest that there are distinct 
emission mechanisms at work and/or that there are 
distinct differences in the plasma conditions in the source 
region. 

Recently, Wong and Goldstein [1990] suggested that 
the bursty radio emissions from planetary magneto- 
spheres could be due to a temperature anisotropic beam 
instability (hereafter called TABI). This mechanism is 
essentially the same as the early proposal by Melrose 
[1976] for the generation of auroral kilometric radiation 
(AKR) and Jupiter’s decametric radiation (DAM), which 
have many similarities to the smooth emissions seen 
at Uranus and Neptune. However, general interest in 
this instability rapidly faded when Wu and Lee [1979] 
developed the theory for the electron cyclotron maser 
instability. This maser instability, which is driven by a 
population inversion in the electron perpendicular velocity 
(i.e., df/dv > 0), is able to provide a ready explanation 
for many of the observed features of AKR and DAM, 
including frequency cutoff of the excited radiation, their 
polarization and brightness temperature, and the overall 
characteristics of the electron distribution in the nightside 
auroral region. 

Nevertheless, there is growing evidence that the types 
of distributions required to drive the TABI can also 
be present, at least on occasions [Menietti and Burch , 
1991]. There is also some recent evidence from laboratory 
experiments that the TABI can be generating bursty 
electromagnetic radiation [Urrutia and Stenzel, 1984; 
Goldman and Newman , 1987], The TABI differs from the 
maser instability in that it is driven by a temperature 
anisotropy where Tj_ Tj| and no explicit positive 
gradients (i.e., inverted populations) need be present. 
It is similar to the maser instability in that it still 
involves generation of electromagnetic radiation via a 
gyroresonant wave-particle interaction. Newman et aL 
[1988] subsequently performed a comprehensive study of 
the mechanism including numerical simulations for the 
case when the ratio of the plasma frequency (wpe) to 
gyrofrequency (O e ) was greater than 2.5. TTon^ and 
Goldstein [1990] have extended the theory to the region 
of ufpe/fle ~ 1 typical of planetary magnetospheres. They 
showed that the TABI can easily generate radiation as 
much as a factor of 2 above the electron cyclotron 
frequency, with the radiation being beamed into a filled 
emission cone of half-angle less than 30° . 

In this paper one-dimensional relativistic particle-in-cell 
simulations are used to investigate the radiation from 
electron beams with a large temperature anisotropy for 
conditions relevant to planetary magnetospheres. The 
use of simulations is critical to the problem because 
these types of distribution are not only unstable to the 
TABI but also to electrostatic beam instabilities which 


generate Langmuir waves. These two instabilities compete 
for and modify the free energy in the distribution as 
they try and reduce the gradients in the distribution. 
This competition between instabilities is neglected in 
all the above theoretical treatments of the TABI. It is 
shown that these beams, which would typically originate 
from the sporadic or impulsive injection of energetic 
electrons, can generate electromagnetic radiation which 
should be able to escape into the solar wind despite 
the growth of the electrostatic instability. This escaping 
radiation is most likely generated under conditions when 
the ratio of the electron plasma frequency u j pe to the 
electron cyclotron frequency is between about 0.3 
and 1.0. The efficiency increases with energy of the 
beam electrons. This parameter regime complements that 
required for the maser instability, which tends to be 
restricted to Wpe/fte < S 0.3. Thus the TABI provides a 
natural explanation for the apparent differences in source 
location and spacecraft event times between the smooth 
and bursty emissions. 

2. Comparison Between Maser 
and Beam Instabilities 

The processes that may be involved in the generation 
of both the smooth and bursty emissions are illustrated 
in Figure 1. The standard model for producing the 
smooth emissions is shown in Figure la. Because of the 
long duration of the smooth emissions, there must be a 
continual resupply of the energetic electrons producing 
the radiation; otherwise the driving mechanism would 
quickly saturate and the emissions would cease. For such 
a continual flow of energetic electrons down the field lines 
a loss cone distribution would form, as electrons with 
large pitch angles mirror and move back up the field lines 
while electrons with small pitch angles precipitate out 
of the flux tube and are lost to the lower atmosphere. 
Such loss cone distributions are unstable to the maser 
instability, which generates the observed radiation for the 
smooth emissions. 

For the bursty emissions, the physics involved is 
substantially different. In particular, the presence of short 
bursts implies that the energetic electrons responsible for 
the emissions are being injected intermittently. For rapid 
or impulsive injection, only certain particles can reach 
a point further down the field line at any particular 
time, since slow particles would not have had sufficient 
time to reach it and fast particles would have already 
passed by [e.g., White et al., 1986; Winglee and Dulky 
1986]. This velocity filtering occurs irrespective of the 
initial distribution, which could be beamlike or just a 
hot Maxwellian of solar wind origin. As a result of 
these time-of-flight effects, the local distribution can be 
narrowly confined in the parallel velocity, allowing the 
distribution to develop a beamlike feature. Furthermore, 
since the time-of-flight effects depend only weakly on the 
perpendicular velocity (assuming only weak magnetic field 
convergence), the particles arriving at any particular point 
can have a much larger spread in their perpendicular 
velocity than in their parallel velocity, effectively giving 
the beam a temperature anisotropy. 

Note that these types of distributions can develop from 
the same initial distribution as that driving the loss cone 
formation; the main difference in the evolution of the 
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Fig. 1. Schematic showing (a) the development of a 

loss cone distribution during continuous injection of energetic 
electrons and the subsequent generation of maser emission 
responsible for the smooth emissions and (b) the development 
of a temperature anisotropic beam during impulsive injection. 
This type of beam can drive a modified Weibel instability 
or TABI which can generate radiation which can possibly 
account for the bursty emissions. 


distribution arises from the injection period relative to 
the transit time of the particles along the flux tube. For 
loss cones to develop, the injected period must be much 
longer than the transit time, i.e., at least greater than 
about a few minutes. Temperature anisotropic beams 
would probably be generated if the injection occurs 
sporadically with each burst having a period smaller 
than about a minute or so. 

When a temperature anisotropic beam develops, a 
modified kinetic “Weibel instability” [Newman et al ., 
1988] or TABI can develop which grows off the difference 
in gradients between the perpendicular and parallel 
velocities. The instability does not explicitly need a 
positive gradient in the distribution (unlike the maser 
instability) although a positive gradient associated with 
the beam component is often assumed in theoretical 
treatments. When the ratio of the plasma frequency 
u/pe to the electron cyclotron frequency 0 C is much 
greater than 1, modified whistlers are excited [Newman 
et a/., 1988], However, if u> pe /Cle < 1, then Wong and 
Goldstein [1990] predict that the dominant emission is 
in the right-hand circularly polarised (RCP) x mode 
with v > O e , and it is this radiation which they 


propose as the possible source of the bursty planetary 
emissions. In the following, the nonlinear evolution of 
the TABI is examined through particle simulations in 
order to compare its characteristics with that of the 
maser instability and, thereby, provide better insight into 
the possible origins of the bursty and smooth emissions. 

3. Particle Distributions and Particle Heating 

The nonlinear evolution of the temperature anisotropic 
beam instability is investigated via one-dimensional (three 
velocity) relativistic electromagnetic particle simulations 
[cf. Newman et a/., 1988]. This code solves self- 
consistently the full set of Maxwell’s equations and allows 
the evaluation of the competition between electromagnetic 
and electrostatic beam instabilities. The system is 
assumed to be periodic with the wave vector h directed 
along the magnetic field, which is in the x direction. 
The size of the system L is assumed to be 1024 Debye 
lengths, and the corresponding wavenumber of a mode 
m in the system is then 2srn/£. 

In all the following simulations the density of the 
plasma is kept constant and the magnetic field strength 
is varied so that v pe / O e ranges from a maximum value 
of 2.5 down to a minimum value of 0.25. This parameter 
regime extends that considered by Newman et al. [1988] 
where only values of w p€ / O e >2.5 were considered. 

The plasma is assumed to consist of (1) a cold ambient 
plasma with a Maxwellian distribution with a thermal 
speed vxc equal to 0.04c, which comprises 90% of the 
total electron density, and (2) a beam component with 
a parallel temperature the same as the cold component 

Early Evolution of the Distribution 


(a) <o pe t = o 


(b) (0 p# t = 60 




Final Distributions 

* (C) <o pe / = 2.5 (d) (0 pe / £l e = 0.4 




“ ° v /v u 

v ll /v Tc 

Fig. 2. The evolution of the total electron distribution 
for v* = 4u Te . (a) The initial distribution, (b) The early 
development of a weak electrostatic beam instability causes 
some parallel heating, (c) and (d) The strongest heating is 
driven by the TABI which produces scattering in both pitch 
angle and energy, as indicated by the arrows. This latter 
heating is dependent on u; M /n« while the heating by the 
electrostatic instability is not. 
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and a perpendicular temperature 36 times the parallel 
temperature, i.e., v^\\e — v Ta v TLE = 6v T||B* Thk ktter 
component makes up the remaining electron plasma 
density, i.e., 10% of the total electron density. These 
parameters are •iTnilnr to those used by Newman et cd. 
[1988]. Three different beam speeds are considered: 
V}, equal to 4, 6, and 8 tj^ c . 

The initial distribution comprising these two compo- 
nents for Vb = 4t >Te ** shown in Figure 2a. It is unstable 
to both the electromagnetic beam instability and a 
weak electrostatic beam instability (similar to the nor- 
mal bump-in- tail instability). The electrostatic instability 
quickly saturates (by Wpet — 60) and produces only some 
flattening of the distribution in the range 3 ~ vu/vx c ~ 4 
(Figure 2b). Most of the particle heating is driven by 
the electromagnetic beam instability, as illustrated in 
Figures 2c and 2d which show the velocity distributions 
at saturation (w pe * ~ 1000) for w pt / 0« equal to 2.5 and 
0.4, respectively. 

The electromagnetic instability acts on the distribution 
in an effort to reduce the temperature anisotropy. For 
utpe/Oe ~ 1 (Figure 2c), this isotropization is achieved in 
two ways: (1) low pitch angle electrons are scattered 
forward to higher t?|| and (2) high pitch angle electrons 
are scattered to lower pitch angles and lower energies. For 
w p «/Oe $ 1, the beam particles primarily experience only 
scattering to lower pitch angles, and the acceleration 
of low pitch angle particles to higher energies is 
significantly reduced. This change in the heating of 
the particle distribution is due to the development of 
different unstable modes as Wpe/fU is varied (see section 

4 >- 

Figure 3 shows the temporal evolution of the distri- 


bution when the beam speed is increased to 6o^ e and 
Wp«/0« is held fixed at 2/3. Because of the higher 
beam speed, the electrostatic instability is very much 
stronger and quickly flattens the beam feature in vp 
as seen in Figures 3a and 3b. As a result of the 
growth of this instability, the distribution no longer 
has a positive slope anywhere. Despite this flattening 
of the beam feature, there are still large temperature 
anisotropies (i.e., df/dv± C df/dv || ) on both the leading 
and trailing edges of the beam feature. It is this 
temperature anisotropy which drives the TABI at later 
times and produces the continued heating of the particle 
distribution (Figures 3c and 3d), similar to Figure 2. 
Note that, unlike the maser instability, a population 
inversion is not required for the growth of the TABI. 

The time history of the energy going into the 
electrostatic waves is shown in Figure 4 for the three 
different assumed beam speeds. The electrostatic wave 
energy for a given beam speed did not vary significantly 
for the different magnetic field strengths (i.e., Wpe/He) 
considered, since the electrostatic waves are generated 
by a Landau resonance and not by a gyroresonance. 
It is seen that the peak energy going into these 
waves increases with increasing beam speed. However, 
in all cases this wave energy is mostly reabsorbed as 
the particle distribution evolves from one which has 
a positive slope in Ojj to one in which such positive 
slopes are absent. The growth and reabsorption of these 
waves occurs on relatively fast time scales, within about 
wpet £ 200. 

The corresponding amount of energy being converted 
to electromagnetic wave energy is shown in Figure 5. The 
electromagnetic wave energy differs from the electrostatic 


^po^e = 0.667, v b /v Tc = 6 
(a) O pe t = 0 (b) <0 pe t = 60 



Fig. 3. The temporal evolution of the electron distribution for v* = 6v Tc and = 2/3. Because 

of the increased beam speed there is much strong parallel heating of the electron distribution via the 
electrostatic beam instability. However, at late times (in Figures 3c and 3d) the TABI is still able to grow 
and produce further particle scattering in pitch angle and energy. 
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Fig. 4. The time history of the particle energy being 
converted into electrostatic wave energy. Irrespective of the 
beam speed, the energy rapidly increases and is then 
reabsorbed. 


wave energy in three ways. First, it grows on much longer 
time scales than the electrostatic wave energy. Indeed, 
there is a tendency for the electromagnetic waves to 
have to wait until the electrostatic waves have saturated 
and been reabsorbed before the TABI becomes intense. 
This delay arises because the fast growing electrostatic 
instability modifies the electron distribution, which in 
turn tends to destroy the slower growing gyroresonant 
interactions. Second, the bulk of the wave energy is not 
reabsorbed. Third, the efficiency depends strongly on the 
value of Wpe/Cle. 

For the lowest beam speed, the efficiency increases with 
increasing u> pe /Q c (Figure 5a). As the beam speed is 
increased to Vf, = 6t>r c (Figure 5b), the overall efficiency 
is larger with peak efficiency increasing about 10% for 
^ 2/3. For Vb = Svj» e (Figure 5c), the marimnm 
efficiency is about the same but occurs at a lower value 
of Wpe/fto equal to about 1/2. This change in the 
efficiency is consistent with the increase in the predicted 
growth rate for increasing beam speed [cf. Wong and 
Goldstein, 1990]. 

The efficiencies shown in Figure 5 of a few percent 
are comparable to that of the maser instability [e.g., 
Pritchett , 1986], In other words, the TABI under certain 
conditions can produce radiation as intense as the 
maser instability. However, the two instabilities occur in 
different but complementary plasma regimes. Specifically, 
the maser instability is know to be relatively strong 
for Upe/Oe ~ 0.3 while the above results show that the 
TABI is primarily effective for Wp e /n e £ 0.3. 

It should be noted that the above simulations treat 
the system as homogeneous so that the true efficiency 
of the instability may be strongly modified by inho- 
mogeneities in the actual source region of a planetary 
magnetosphere. This is also true for the maser instabil- 
ity. For example, wave propagation out of the source 
may limit the amplification of the wave. On the other 


hand, particle propagation effects may help maintain the 
ini tial distribution, as heated particles move out and new 
particles propagate into the source region. As a result, 
the instability may be sustained over many e- folding 
periods rather than quickly reaching saturation as the 
gradients in the distribution are eroded. Evaluation of 
these competing effects needs detailed information about 


the source region which is currently not available. 

4. Characteristics op thb Wave Spectrum 
The characteristics of the excited waves are expected 
to change as Wpe/fte is varied from above to below unity. 




CO 


pe 


t 


Fig. 5. The time history of the energy going into the 
electromagnetic wave for the different values of beam speed 
and u> 7 «/n«. As the beam speed is increased beyond a few 
thermal speeds, the efficiency increases to a maximum of 
about 10%. This peak efficiency is, however, attained over 
longer periods and lower values of u; p# /n«. 
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In particular, for Wpe/fle ~ 1, the whistler can have a 
frequency comparable to the electron cyclotron frequency 
due to a Doppler shift associated with the presence 
of the beam. This Doppler shift allows gyroresonant 
interactions to occur. However, when utpe/fte £ 1* the 
mn^mnm frequency of the whistler is limited to u/ pe , 
where it becomes increasingly electrostatic so that 
gyroresonant interactions are suppressed. However, the 
cutoff for the x mode and resonance for the s mode 
start to approach the cyclotron frequency, and it is 
these modes which can then become unstable through 
gyroresonant interactions if the beam velocity is large 
enough. This change in the mode structure is illustrated 
in this section. 

Figure 6 shows the spectrum derived from integrating 
the power in the different k modes over the full 
period of the simulations for v^jv^c = 4. In order to 
separate spectra from different modes, some of the lines 
have been dotted. It is seen that, for w P c/Oe = 2.5, 
the dominant emission is at about 0.5 Qe and is 
left-hand polarised (LCP). This mode is essentially a 
beam-modified backward propagating whistler. Because 


its frequency is below He, this emission is expected to 
be trapped in a planetary magnetosphere. 

At the intermediate values of u> pe /Oe (Figures 6b and 
6c), there is still relatively strong emission at 0.5 O e . In 
addition, strong emissions appear at frequencies slightly 
above O tf . This is the RCP x mode predicted by Wong 
and Goldstein [1990], and it is this radiation, particularly 
for the lower values of « pe / O e , which can possibly escape 
into the solar wind and produce the bursty emissions. 
However, in the present case the beam speed is small, 
and most of the wave energy has frequencies below 
the local x mode cutoff and cannot directly escape the 
plasma. Only at the lowest value of w p «/Oe shown is 
there a small tail (ft! 1% of the electromagnetic energy) 
with frequencies above the local x mode cutoff. 

Nevertheless, there is still a possibility that some of 
the radiation can escape because the electron injection 
must be associated with development of field-aligned 
currents. These currents can produce perturbations in 
the magnetic field of a few percent of the local ambient 
magnetic field in the nightside auroral region, and there 
is no reason to believe that such magnetic perturbations 





(O/flg 

Fig. 6. The frequency spectrum of the different fc modes when v./vt. = 4 oud for the four different 
Some of the spectra appear as dotted lines in order to highlight differences between modes. For 
£ 1, the do minan t mode is LCP. However, as is decreased, the RCP waves become relatively 

more intense with a significant portion above O. but below the x mode cutoff. Only in Figure fid is there 
a weak tail in the spectrum with frequencies above the x mode cutoff. 
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should not also be present in the magnetospheres of the 
outer planets. Such perturbations, could modify the local 
cyclotron and cutoff frequencies so that they temporarily 
fall below the frequency of the excited radiation. In this 
case, inode conversion of the induced TABI radiation 
can allow the escape of some of the radiation, even for 
this relatively low beam speed. 

As the beam speed is increased, a huger portion of 
x mode waves is generated with frequencies above the 
local x mode cutoff. As an example, Figure 7 shows the 
spectra for v^/vxt = 6 and for three values of uip«/ 0« 
where the relative efficiency of the TABI is near its 
maximum. It is seen that emissions in the LCP mode are 
all but suppressed and that the dominant emissions are 
in the RCP mode with frequencies above the cyclotron 
frequency. Moreover, for the smallest value of Upe/O e 
shown, the dominant emission has a frequency above the 
x mode cutoff which is denoted by ui® in the figure. 

Further increases in the beam speed (Figure 8) produce 



CO/fl, 


Fig. 7. As in Figure 6 but for vt/v T « = 6. For this higher 
beam speed the radiation with frequencies above the x mode 
cutoff hi. actually dominates the spectrum at the lowest value 
of ai,«/Q« shown. 



co/Q* 


Fig. 8. As in Figure 6 but for v»/vt c = 8. There is a 
further upshift in the frequency spectrum with more of the 
radiation appearing above oi., and the growth of all the LCP 
waves is suppressed. 

a continuing upshift in the frequency of the dominant 
modes. For the highest value of co pe /n e shown (Figure 
8a), the upper cutoff of the excited radiation is nearly 
1.5 O e compared with only about 1.17 H e in Figure 6. 
At smaller values of w ?9 f{\ e a growing fraction of the 
wave energy is above the x mode cutoff. 

Note that while the frequency of the excited radiation 
is increasing with beam speed, it increases at a smaller 
rate than predicted by Wong and Goldstein [1990]. This 
difference arises from the presence of the fast growing 
electrostatic waves which are not included in the linear 
theory of Wong and Goldstein [1990]. These waves tend 
to cause a reduction in the average drift speed of the 
energy beam particles and hence in the Doppler shift. 

5. Conclusions 

One-dimensional electromagnetic particle simulations 
have been used to evaluate the characteristics of radiation 
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from electron beams with temperature anisotropies in 
order to identify the origin of bursty planetary emissions. 
Such beams can arise for example during the impulsive 
or sporadic injection of energetic electrons in which 
time-of-flight effects restrict the parallel velocity of the 
particles arriving at any specific point and time down 
the field lines. Terrestrial observations show that such 
temperature anisotropies in association with electron 
beams can indeed be present at least in the Earth’s 
auroral region [cf. Menittti and Burch, 1991]. In 
contrast, the smooth planetary emissions are believed to 
be generated by the cyclotron maser instability (CMI), 
which is presumably driven by the continuous injection of 
the energetic electrons and the subsequent development 
of loss cone distributions. 

These temperature anisotropic beams are unstable to 
a modified kinetic Weibel or temperature anisotropic 
beam instability (TABI). For w pe / 0 e Z 1, most of the 
wave energy generated is LCP with frequencies below the 
cyclotron frequency and is therefore trapped. However, 
for Wpe/fte £ 1 the dominant mode is RCP, either in 
the form of a modified whistler or in the x mode. The 
latter radiation is strongest for 0.25 £ Upe/H* ~ 0.5 and 
where beam speed is much greater than the ambient 
thermal velocity. Under these conditions, the excited 
radiation can directly escape from the plasma and be 
detected remotely by a spacecraft. At moderate values 


of Upe/fle (i.e., between about 0.5 and 1.0), much of the 
excited radiation is below the x mode cutoff, and the 
only possibility for the radiation to escape is via mode 
conversion across adjacent regions where the magnetic 
field is locally lower. Such perturbations in the magnetic 
field may indeed be present since field-aligned currents 
are likely to be associated with the electron injection. 
Detailed modeling of these latter processes still needs to 
be performed in order to establish the exact maximum 
value of Upe/n e for which the TABI radiation will be 
trapped in planetary magnetospheres. 

In addition, the TABI is able to convert particle energy 
into electromagnetic energy with similar efficiencies as the 
maser instability. As a result, the escaping radiation from 
the TABI, particularly when wpe/He — 0.4, should have 
a comparable intensity to that of the maser radiation, 
which is strongest in regions where w pc /Oe £ 0.3. 

Thus the above features of the TABI radiation 
strongly complement those of the maser radiation, both 
in the plasma conditions in which escaping radiation 
is generated and in the injection of the energetic 
electrons as it changes from continuous to sporadic. This 
complementary nature provides a ready explanation for 
the apparent offsets in source location and spacecraft 
event time of the smooth and bursty emissions, as well 
as their different spectral characteristics. 

The model presented here suggests a radiation pattern 


Clfl (SMOOTH) 
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and source location similar to those shown in Figure 
9, with the source regions of the two emissions ap- 
proximately adjacent , as Upe/He and the characteristics 
of the electron injection vary across the polar regions. 
As a result of the different emission mechanisms and 
their source location, a spacecraft would be expected 
to see the two sides of the hollow cone of the smooth 
emissions followed by a single patch of bursty emission 
associated with emission from a filled cone, similar to 
the observations of Farrell and Calvert [1989]. 
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Abstract We show that an energetic electron distribution which has a temperature 
anisotropy ( T ^ > 7^), or which is gyrating about a DC magnetic field, can generate 
electron cyclotron waves with frequencies below the electron cyclotron frequency. 
Relativistic effects are included in solving the dispersion equation and are shown to be 
quantitatively important. The basic idea of the mechanism is the coupling of the beam 
mode to slow waves. The unstable electron cyclotron waves are predominantly 
electromagnetic and right-hand polarized. For a low-density plasma in which the electron 
plasma frequency is less than the electron cyclotron frequency, the excited waves can have 
frequencies above or below the electron plasma frequency, depending upon the parameters 
of the energetic electron distribution. This instability may account for observed Z mode 
waves in the polar magnetosphere of the Earth and other planets. 


Introduction 

The Earth’s polar magnetosphere has long been recognized 
as an active region of wave activities. In the last two decades, 
numerous spacecraft have sampled this region of space and 
have identified a large variety of wave modes, most notably the 
auroral hiss and the auroral kilometric radiation (AKR). For a 
review, see Shawhan [1979]. These waves are believed to play 
an important role in various plasma processes in the magneto- 
sphere through wave-particle interactions. Examples of such 
processes include the diffusion of auroral electrons by electro- 
static electron cyclotron waves [Kennel and Ashour-Abdalla , 
1982], the generation of AKR through relativistic cyclotron 
resonance [ Gumett , 1974; Wu and Lee , 1979] and the accelera- 
tion and heating of ions and electrons by waves in the auroral 
region leading to the formation of ion and electron conical dis- 
tributions (see, for example, Chang et ai [1986], Crew et ai 
[1990], Lysak [1986], Temerin and Cravens [1990], and Wong 
et ai [1988], among others). 

The whistler mode has been studied extensively in the past, 
primarily due to the interest in VLF emissions, auroral hiss, 
and lightning-related phenomena. In the midaltitude polar 
magnetosphere, the electron plasma frequency 
G) e =(4m e e 2 fm e ) 112 is typically less than the electron 
cyclotron frequency Q e = \e\B 0 /m e c f where m e is the electron 
rest mass, n e is the electron density, B 0 is the magnitude of the 
ambient magnetic field, and c is the speed of light. Under this 
condition, the whistler mode propagates between the lower 
hybrid frequency, which is approximately the ion plasma fre- 
quency CO,, and co*. The frequently observed whistler mode auro- 
ral hiss is believed to be generated near the resonance cone 
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either by precipitating electrons or by upward moving electron 
beams [Gumett et ai , 1983]. The whistler waves excited near 
the resonance cone are quasi-electrostatic waves with negligi- 
ble magnetic components; however, analysis of the “funnel- 
shaped” auroral hiss observed by DE 1 indicates that the auroral 
hiss has a considerable magnetic component [Gumett et ai , 
1983], Subsequent stability analysis using the observed parti- 
cle distributions has revealed that the electron acoustic wave, 
rather than the whistler wave, is the dominant unstable wave 
driven by such electron beam distributions [Lin et al. , 1984, 
1985; Tokar and Gary , 1984]. The electron acoustic wave is 
electrostatic in nature, and thus might account for the electro- 
static component of auroral hiss; however, the origin of the 
electromagnetic component of the hiss still remains unan- 
swered. 

In a different context, Benson et al [1988] have reported 
ground-based detection of waves in the frequency range 150— 
300 kHz, indicating the generation of field-aligned waves in 
the auroral zone. The field-aligned waves observed by Benson 
et al [1988] are electromagnetic and fall into the frequency 
range of the whistler mode and overlap the frequency range of 
AKR. However, because the theories usually proposed for the 
generation of AKR produce radiation in the extraordinary and 
ordinary modes which cannot reach the ground, the field- 
aligned radiation observed by Benson et al [1988] must be 
generated by a different mechanism. Motivated by these 
observations, Wu et al [1989] showed that an energetic elec- 
tron population with either a temperature anisotropy or a 
trapped type distribution can generate field-aligned waves in 
the observed frequency range. More recently, Ziebell et al 
[1991] have generalized Wu et al's [1989] theory to investi- 
gate propagation effects on the amplification of radiation 
along the field line. 

Previously we have studied the generation of bursty radio 
emission by low-density anisotropic or gyrating electron 
beams [Wong and Goldstein , 1990] (hereinafter paper 1). In 
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the present study we extend that analysis by including rela- 
tivistic effects and show that the frequency range of the unsta- 
ble waves excited by dense gyrating and/or anisotropic distri- 
butions can extend below ft, or even below to,. In this regard 
our analysis is a generalization of the previous work by Wu et 
al. [1989] because we extend into the relativistic regime, 
include oblique propagation, and use a more general distribu- 
tion function that includes a beam component. The process we 
describe appears to provide an alternative to the cyclotron 
maser instability for generating Z and O mode radiation. 


Physical Model and Results 

The physical model we consider is the same as previously 
discussed in paper 1: a gyrating or anisotropic electron beam 
in a cold background magnetized plasma. Because the frequen- 
cies we consider are close to ft, or to,, we can ignore the ion 
contribution to the plasma dielectric. In paper 1 we discussed 
the conditions under which such distributions can arise and we 
do not repeat that discussion here. The general form of the 
energetic electron distribution we use is 


A( w lh w l)“ 3/2 2 eX P z 

k a jjy ajjj 


(“L-“10) 2 («ll~“fe) 2 


«ll b 


( 1 ) 


where u = p /m, is the momentum per unit mass, = 

(2icT 1 ^//n <r ) 1/2 , K is Boltzmann’s constant, u b and u ±0 are the 
momenta per unit mass of the electron beam parallel and per- 
pendicular to the magnetic field, respectively, and the normal- 
ization constant is found from A -1 = exp(-M X ^/a 1 ^ 2 )[l-i(M JL0 
/ajj, ) Z(i(u 10 /an , )] to ensure that j (u) = 1 . 

A general derivation of the relativistic electromagnetic dis- 
persion equation for arbitrary directions of propagation can be 
found in the work by Baldwin et al. [1969]. In the appendix we 
give a brief description of how we solve the dispersion equa- 
tion numerically. Most of the discussion below will be limited 
to parallel propagating solutions, primarily because the maxi- 
mum growth is generally at 0°. However, we have also solved 
the more general obliquely propagating relativistic dispersion 
equation, and we discuss below some of those results (Figure 6). 

For comparison purposes, we first show solutions to the 
nonrelativistic dispersion equation for parallel propagating 
modes excited by trapped electrons with energies of the order of 
several keV, a case discussed previously by Wu et al. [1989]. 
The dispersion relation was solved using n b = 0.9, n e = 9n c , 
where n b and n c are the densities of the energetic and back- 
ground electrons, respectively; and n e = n b + n c . Note that, in 
contrast to paper 1, we emphasize here a regime in which the 
density of the energetic electron component, n b , exceeds that 
of the background, n c . We used five values of co,/ft,, viz. 0. 1 , 
0.2, 0.4, 0.6, and 0.8. In the example shown in Figure 1, the 
drift of the electrons is zero, and the parallel and perpendicular 
temperatures of the energetic electrons are both 2.5 keV. The 
ring speed is m X q = 2.0 in units of a j^. In Figure 1 we plot the 
real and imaginary frequencies (co r and to,*, respectively) as func- 
tions of wave number. The frequencies are normalized by ft,, 
and the wave numbers by ft,/c. Note that to r and oo,- are plotted 
on linear and logarithmic scales, respectively. One should 
keep in mind in these analyses that because the growth rates are 
functions of co,/ft, and n b /n e , specific values of density need 
not be specified. 

As is evident from Figure 1, the instability exists for the 
entire range of the values of to,/ft, considered, as does the band- 


width of the instability. For small values of (0 ,/ft,, the fre- 
quency of the unstable waves lies between to, and ft,, which is 
in the frequency regime of the Z and O modes. For larger values 
of to,/ft„ the frequency of the unstable waves can extend below 
to,; the waves are then in the whistler and Z mode branches for 
frequencies above the left-hand cutoff. This result is essen- 
tially the same as was found by Wu et al. [1989], as can be veri- 
fied by comparing our Figure 1 with their Figure 3. 

In Figure 2 we show the effect on the magnitude of the growth 
rates when using the relativistic dispersion relation. The 
parameters are the same as we used in Figure 1, except that now 
co ,/ft, = 0.8. It is clear that both the bandwidth and magnitude 
of the instability are reduced. The reduction in bandwidth is 
from A£c/ft, = 0.6 to Afcc/ft, = 0.4, while the maximum 
growth rate drops by nearly a factor of 5. Thus it is clear that 
relativistic effects can be quantitatively important in this 
class of instabilities. For the remainder of this paper, all 
results will be computed using the relativistic dispersion equa- 
tion. 


Parallel Propagating Electron Cyclotron Wave 



Figure 1. The unstable roots of the nonrelativistic 
dispersion equation plotted for comparison with Figure 3 of Wu 
et al. [1989]. The distribution function of the energetic 
electrons, given by (1), is normalized so that n b = 0.9 and n e = 
9 n c . The parallel and perpendicular temperatures of the 
energetic electrons both equal 2.5 keV, and the ring speed is 
“LcJQ-Jdb ~ 2-0 (the beam speed u\\ = 0). The real and imaginary 
parts of 0) are plotted on linear and logarithmic scales, 
respectively, for ( 0 ,/ft, = 0.1, 0.2, 0.4, 0.6, and 0.8. 
Frequencies are normalized by ft, and wave numbers by Q.Jc. 


Parallel Propagating Electron Cyclotron Wave 



Figure 2. Comparison between solutions of the relativistic 
and nonrelativistic dispersion equations for the same 
parameters are were used in Figure 1 except that only the results 
for co,/ft, = 0.8 are shown. 
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For zero drift speed (wj_o = u b = 0), this instability can also be 
driven by a temperature anisotropy, as is illustrated in Figure 3. 
Here the density ratio between the energetic and background 
electron distributions remains as before, ( o<JQ e = 0.6, and the 
parallel temperature of the energetic electrons is kept constant 
at 2.5 keV, while the perpendicular temperature is increased 
from Tjjy/Tfti, = 4 to T Lb /T^ = 10. For Ti b /T\$ = 4 the growth 
rate is very small (co, mflJt /D, = 10~ 4 ), but it increases by more 
than 2 orders of magnitude as the anisotropy increases to 
T ±j>/T\\b = Clearly, at T lb /T^ = 4 the instability is near 
threshold (note the rapid increase from T^/T^ = 4 to T Lb /T^ = 
5) — subsequent increases in T produce far smaller 
changes in the maximum rate of growth. 

Figure 4 illustrates the behavior of the instability as the ring 
speed is varied from = 2.8 to 3.2. The parameters used 

in these calculations are again the same as above, except that 
w e /£l e = 0.4. The growth rate is a fairly sensitive function of 
w ±c/ a ±i>’ which can be understood because the free energy of the 
ring is equivalent to a temperature anisotropy, as discussed by 
Wong and Goldstein [1987], Recall from Figure 2 that rela- 
tivistic effects can significantly reduce the growth rate from the 
values calculated nonrelativistically, especially for relatively 
small values of and/or w 10 . This must be kept in mind 

when comparing Figures 2 and 4 with the nonrelativistic calcu- 
lation shown in Figure 1. 

In Figure 5 we investigate how finite beam speed affects the 
instability as u b fa ^ is varied from 0 to 2. The parameters used 
are again as in Figure 1 , except that = 2.5 and = 


0.6. Small increases in u b /c^ b (from 0 to 1) cause the instabil- 
ity to broaden in bandwidth and to increase slightly in maxi- 
mum growth. Further increase in u b la^ to 2 produces a reduc- 
tion in growth and in bandwidth from what it was at u b /o = 1; 
however, it is still significantly broader than it was initially 
when u b /CL\\ b =0. This behavior contrasts sharply from the 
beam-driven radio emission discussed in paper 1 (cf. Figure 2 in 
that paper) in which the growth rate increased monotonically 
with increasing beam speed. 

We have also investigated the behavior of this instability 
for oblique propagation. Again we take the thermal energy of 
the energetic electrons to be 2.5 keV, the density ratio of the 
energetic and cold components of the electron distribution 
equal to 9, the ring speed = 2.5, and oyQ, = 0.6. The 

distribution has no drift ( u b = 0). The real and imaginary parts 
of co are plotted in Figure 6 as a function of 0, the angle 
between k and B. The plot is constructed at the wave number of 
maximum growth of the instability at 0 = 0°, kcJ& e = 1.28. 
For this wave number, the instability extends to approximately 
10°, and then drops rapidly to zero. This does not necessarily 
mean, however, that there is no instability at larger angles. 
For the parameters of this example, as the propagation angle 
changes, the wave number of local maximum growth at that 
value of 0 shifts to lower values of k, and the cone of unstable 
wave vectors extends somewhat beyond 15°. Nonetheless, the 
instability still grows most rapidly at 0 = 0°. There is little 
Landau damping at these frequencies, and these oblique waves 
are primarily electromagnetic in nature. 


Electron Cyclotron Instability — Variation with Temperature Anisotropy 



Figure 3. Variation in the growth rate of the electron 
cyclotron instability with changes in temperature anisotropy. 
The beam speed is 0, and (o e /Q e = 0.6. The temperature 
anisotropy varies from T^/T^ = 4 to 7^/7]^ = 10. At Tj^/T^ = 
4 the instability is near threshold. 


Parallel Propagating Electron Cyclotron Ware (Relativistic) 



Figure 4. Variation in the growth rate of the instability for 
(Og/Qg = 0.4 as the ring speed is varied from u x0 / a lb = 2.8 to 
3.2. All other parameters are the same as were used in 
Figure 1. 


Electron Cyclotron Instability — Variation with Beam Speed 



Figure 5. The effect of finite beam speed on the instability. 
In this case, u b / is varied from 0 to 2; = 2.5, and 

co elQ e = 0.6. All other parameters are the same as were used in 
Figure 1. 


Obliquely Propagating Electron Cyclotron Wave Instability 
(Relativistic) 



Figure 6. The behavior of the instability for off-angle 
propagation. The thermal energy is again 2.5 keV, u Lo l a ±b = 
2.5, and co g/£l e = 0.6. The beam drift speed u^ b = 0, and kcf£l e = 
1.28. The density of the beam is the same as was used in all 
figures above. 
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As the value of (ojilg increases, e.g., to */£2 r = 0.8, the situa- 
tion changes significantly. For large values of oyO,, maxi- 
mum growth occurs at oblique propagation and at larger values 
of k than for the local maximum at 0 = 0°. We suspect that this 
change in the nature of the solution arises from the fact that the 
real part of the wave frequency is above the electron plasma fre- 
quency when the ratio is small, whereas when co g/Q e = 

0.8, the excited wave has a frequency below co tf ( cf. Figure 1). 
Thus at large values of co^/D^, the excited wave can be in the Z 
and whistler modes, while for small values of (o ^£2*, the wave 
can only propagate in the Z and O modes (for a discussion of 
which wave modes exist at particular frequencies in the auroral 
region, see Figure 5 of Gumett et al. [1983]). This behavior of 
the instability with changing <a e /Q e suggests that the trapped 
(or ring) distribution is an effective population for generating 
whistler waves if the physical parameters are appropriate. 

Discussion 

The fundamental result of this calculation is the demonstra- 
tion that a temperature anisotropy or a ring-beam electron dis- 
tribution can generate electron cyclotron waves with frequen- 
cies that extend from below to above co f , but still below Q e for 
an underdense plasma. The calculation is relativistic and 
allows for arbitrary direction of propagation. Previous work 
[Wu et al., 1989] was confined to parallel propagation and was 
a nonrelativistic calculation. Our result for parallel propaga- 
tion agrees qualitatively with Wu et al.’s results; however, as 
shown in Figure 2, the effect of the relativistic corrections to 
the dispersion equation is to reduce substantially both the 
magnitude of the growth rate and the bandwidth of the unstable 
waves. 

We also investigated the effect that Finite beam speed has on 
the instability, an effect not considered previously. For small 
beam speeds, the growth rates can be enhanced and the band- 
widths broadened. But as the beam speed increases, the growth 
rate decreases (cf. Figure 5). In paper 1 we studied the genera- 
tion of bursty radio emission and showed that in plasma 
regimes in which co*/D f < 1, anisotropic electron beams, or 
gyrating electron beams, can excite directly right-hand polar- 
ized broadband electromagnetic radiation in the X mode. The 
present analysis extends that work to waves with frequencies 
below and also includes relativistic effects in the dispersion 
equation. In addition, the mechanism described in paper 1 gen- 
erated unstable X mode waves when the electron distributions 
had large drift speeds. In contrast, for the parameters used here, 
the present mechanism is suppressed for large electron drift 
speeds and excites waves primarily in the Z or O mode. 
Nonetheless, the free energy source in both situations is the 
temperature anisotropy or ring component of the energetic 
electrons (the present mechanism operates even in the absence 
of any drift). 

Recently, Winglee et al. [1992] have performed a one- 
dimensional (with three components of velocity), relativistic 
electromagnetic particle simulation that showed that an 
anisotropic electron distribution can generate freely propagat- 
ing right-hand circularly polarized electromagnetic radiation if 
the beam speed is sufficiently large — in agreement with the 
results of paper 1. Winglee et al.’s simulations are periodic, 
and the wave vectors are parallel to the magnetic field. Besides 
the freely propagating wave, Winglee et al. [1992] also found 
an instability with frequencies below £2 e — a result consistent 
with the theoretical analysis here. For lower values of the 
beam speed they found that the dominant wave mode consisted 
of waves with frequencies co r < whereas for large beam 


speeds, the dominant mode was the freely propagating X mode, 
also in agreement with the results shown here and in paper 1. 

The motivation for the present work was observations of Z 
mode radiation in the Earth’s auroral zone [Gumett et aL , 1983] 
and similar phenomena at other planets [Gumett et al ., 1990; 
Kurth and Gumett , 1991]. Although Z mode radiation can be 
generated by the cyclotron maser instability (see, for example, 
Hewitt et al. [1983]), ray tracing studies [Menietti and Lin , 
1986] indicate that cyclotron maser resonance may not be the 
only mechanism responsible for the observed Z mode emission 
because the generation and subsequent refraction of the waves 
do not appear consistent with observations. Electron 
cyclotron waves generated by temperature anisotropies or ring- 
beam distributions can have frequencies between co tf and in 
the Z and O mode ranges. Because the waves are mostly right- 
hand polarized, we suspect that this mechanism mainly pro- 
duces Z mode radiation. The electron cyclotron waves excited 
in this frequency range are predominantly field-aligned and 
have a much broader bandwidth in frequency than waves gener- 
ated by the cyclotron maser instability. Thus the electron 
cyclotron wave generation mechanism discussed here may pro- 
vide another, complementary, explanation for observed Z mode 
radiation. Because a ring-beam distribution is also capable of 
generating Z mode radiation via the cyclotron maser instabil- 
ity, determining which mechanism actually dominates in any 
specific situation requires detailed modeling that is beyond the 
scope of this paper. 

Conclusions 

We have examined the role of relativistic effects on the gen- 
eration of electron cyclotron radiation by anisotropic energetic 
electrons. The anisotropies are assumed to be due to either tem- 
perature anisotropies or ring (trapped) distributions. Although 
our results are qualitatively in agreement with previous work 
[Wu et al ., 1989], relativistic effects tend to suppress signifi- 
cantly the growth rate. We have also found it possible to 
excite electromagnetic whistler waves with maximum growth at 
a substantial angle from the background field for large values of 
co eJQg (~ 0.8, but still less than unity). The excitation mecha- 
nism discussed in this paper may provide an alternative expla- 
nation to the cyclotron maser instability for the generation of 
Z mode radiation observed in planetary magnetospheres. 

As shown in paper 1 and in the simulations of Winglee et al. 
[1992], a gyrating or anisotropic electron distribution can 
excite a variety of waves in several different frequency ranges. 
In addition, the same class of distributions is also unstable to 
the cyclotron maser instability. Which mechanism (or wave 
mode) dominates will depend on the detailed form of distribu- 
tion and on details of the background plasma parameters. A 
comprehensive comparison of the relative importance of these 
mechanisms will have to await more extensive analytic work 
coupled with two-dimensional numerical simulations, includ- 
ing a study of the convective properties of these instabilities 
along auroral field lines (for example, Ziebell et al. [1991]). 

Appendix 

In this derivation, two electron components are immersed in 
a uniform, external magnetic field B = B 0 e z : a low-energy 
thermal background and an energetic relativistic population. 
The density ratio between these two populations can be arbi- 
trary. Because our interest is confined to waves with frequen- 
cies of the order of the electron cyclotron frequency, which is 
much greater than the ion plasma frequency, ion dynamics can 
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be neglected. Consequently, the dispersion equation of the lin- 
earized Vlasov-Maxwell equations can be written in the form 
[Wong et at,, 1989] 


D(k,<o) = 


k 2 r 2 \ c 2 

1-— T l + £ Tkk + IQ“(tt») = 0 (Al) 


a r 


co A 
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where y= (1 + « 2 /c 2 ) l/2 , and a = b,e refers to the background 
and energetic electrons, respectively. The wave vector k is 
assumed to lie in the jc-z plane so that k = kj L e x + k^e z \ I is the 
unit dyadic; co\ = 4nn a e 2 / m a \ u is the momentum per unit 
mass and is given by u = p lm e ; and n a is the number density of 
species a. The tensor T n is defined by 


T n = 


■ 

ti 

inQ e 


*± 

nCl e 


“ii A 


in£l e 


“i KK 




U±J'n 2 


-i“i\U±J n J' n 


^ u, / 2 

iu\\u t J n J’ n 


(A3) 


where J n =J n (k 1 uJQ e ). 

The method for obtaining numerical solutions to the rela- 
tivistic dispersion equation is cumbersome, and here we only 
outline the approach employed in reaching the solutions given 
above. The relativistic factor y is expanded to Oiu^/c 4 ), except 
for the resonant denominator in which we keep the entire 
expression yto - n£2 e - &||U||. One then integrates over Uj| analyt- 
ically and writes the resulting expression in terms of the 
plasma dispersion function (the Z function defined by Fried and 
Conte [1961]). In contrast to nonrelativistic calculations, the 
argument of the Z function now depends on u ± because of the 
factor of y, which appears in the resonance condition. The 
remaining integration over is performed numerically using a 
standard root-finding algorithm to compute (complex) <o. 
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Abstract . "High-latitude” radio emission was briefly observed as Voyager approached Neptune's north 
magnetic pole. The times of maximum emission and of cutoff depend on frequency. Coupled with a 
magnetic-field model, these data yield information on radio-source location and emission angle. In a 
dipole field, L-shell and magnetic longitude locate candidate sources of each frequency. Emission 
angles are then identified with the angle between the modeled magnetic-field direction at the source and 
the line of sight to Voyager 2 at the times of maximum. At each value of L in the range 6 < L < 9, 
there is one source longitude for which emission angle varies smoothly with frequency. The data thus 
limit source locations to a defined set in the dipole field. Examination of more complex magnetic 
models [Connemey et al., 1991] shows (1) assumptions of the dipole-based analysis are inconsistent 
with complex models; (2) the parametric analysis developed here could under some circumstances be 
extended to a more complex model; (3) the model describes constant-field surfaces with steep slopes 
that explain qualitatively the observed features of the radio cutoff; and (4) radio emission occurs in 
association with relatively small-scale magnetic structure described by multipole terms in the magnetic 
model. 


Dipole-Based Analysis 

At Neptune the Planetary Radio Astronomy (PRA) experiment observed smooth radio emission in the 
range 20 to 1326 kHz. Presumably this is extraordinary-mode (X-mode) radiation at the electron 
cyclotron frequency, f c [kHz] = 0.028 B[nT], corresponding to magnetic field B of 714 to 47,357 nT. 
Radio sources at each frequency lie on a surface of constant field. The magnetometer on Voyager 
measured fields as strong as 9,695 nT (271 kHz). Voyager must have passed through each source surface 
with f c < 271 kHz. At each frequency in the range 39.6 to 462.0 kHz reversal of sense of circular 
polarization accompanied an abrupt drop of intensity. This event was interpreted as X-mode occultation 
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Abstract We describe procedures for analyzing data from the Planetary Radio Astronomy 
experiment. A method to deduce the location and characteristics of the radio sources of Neptune's 
smooth recurrent radio emission minimizes limiting assumptions and maximizes use of the data, 
including quantitative measurement of circular polarization. A preliminary broad survey examines 
all possible source locations at key times. Detailed study of specific sources simulates intensity 
and apparent polarization of their integrated emission for an extended time period. Time series are 
modeled here for broad and beamed emission patterns, at two frequencies that exhibit different 
time variation of polarizadon. These incomplete results are based on a dipole magnetic model and 
include only qualitative consideration of polarization transfer. Consideration of a more complex 
model of Neptune's magnetic field overturns basic assumptions of the dipole-based analysis . Both 
magnetic polarities occur in both magnetic hemispheres so either sense of polarization can arise in 
either hemisphere. Irregular elevation of the source surface modifies the radio horizon. Estimating 
the radio effects, we note evidence of association of sources of Neptune’s smooth recurrent 
emission with magnetic multipolar structure. Multipoles and smooth recurrent emission are 
simultaneously visible from Voyager, multipolar fields can account for the observed range of radio 
frequency and for the observed sense of circular polarization. 

1. Introduction 

With a model of the magnetic field, frequency defines the radial coordinate of a radio source. A cutoff 
occurs where the index of refraction vanishes: (f p If) 2 + f c /f = 1, where f c is the local value of the 
electron cyclotron frequency and f p is the plasma frequency. When f p « f c , as required by the 
cyclotron maser theory, the wave propagates at a frequency just above f c (kHz) = 2800 B (gauss). Radio 
sources of a given frequency lie on the corresponding constant-field surface defined by the magnetic 
model. 

The Planetary Radio Astronomy (PRA) experiment measured the time variation of intensity and 
circular polarization over a range of radio frequencies at Neptune [Warwick et al . , 1989]. PRA’s 
incomplete set of measurements precludes direct determination of source direction and polarization at a 
given moment. PRA experimenters have devised a number of ways, listed below, to locate sources and 
estimate the emission angle using the existing data. 

Source location basically assumes extraordinary-mode (X-mode) emission and matches magnetic-field 
direction to the observed sense of circular polarization. Planetary radio sources have traditionally been 
located in the appropriate auroral zone, with emission in a hollow cone at wide angle to the local magnetic 
(dipole) field. Another basic technique identifies the beginning and end of a period of emission with 
source passage across the radio horizon, the limit of visibility of X-mode emission [Desch and Kaiser , 
1987]. A third method depends on the fact that circular polarization detected by the PRA instrument 
vanishes when the source is on the antenna electric plane [Warwick et al ., 1977]; a spacecraft maneuver 


